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We study the electron transport through a system of two 
low-capacitance metal islands connected in series between two 
electrodes. The work is motivated in part by experiments 
on semiconducting double-dots, which show intriguing effects 
arising from coherent tunneling of electrons and mixing of 
the single-electron states across tunneling barriers. In this 
article, we show how coherent tunneling affects metallic sys- 
tems and leads to a mixing of the macroscopic charge states 
across the barriers. We apply a recently formulated RG ap- 
proach to examine the linear response of the system with high 
tunnel conductances (up to 8e 2 /h). In addition we calculate 
the (second order) cotunneling contributions to the non-linear 
conductance. Our main results are that the peaks in the linear 
and nonlinear conductance as a function of the gate voltage 
are reduced and broadened in an asymmetric way, as well as 
shifted in their positions. In the limit where the two islands 
are coupled weakly to the electrodes, we compare to theoret- 
ical results obtained by Golden and Halperin and Matveev et 
ai. In the opposite case when the two islands are coupled 
more strongly to the leads than to each other, the peaks are 
found to shift, in qualitative agreement with the recent pre- 
diction of Andrei et al. for a similar double-dot system which 
exhibits a phase transition. 
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I. INTRODUCTION 

Electron transport through small metal islands 
displays single-electron effects such as Coulomb- 
blockade and gate-voltage dependent oscillations of the 
conductanceETEl. Most of these effects can be explained in 
terms of lowest-order perturbation theorycl, but recently 
it has been found - in experiment and theory - that in 
several regimes higher-order tunneling processes need to 
be taken into account. These are the Coulomb- blockade 
regime, where sequential tunneling is exponentially su- 
pressed, and the regime close to resonances. In the for- 
mer the dominant contribution to the current is due to 
cotunneling processes, where electrons tunnel through 
the system via virtual intermediate states3'u. Close to 
the resonances, cotunneling and higher-order processes 
can significantly affect the linear and nonlinear conduc- 
tance even -though the lowest-order processes are not 
supressedJQ. 

In this work we study equilibrium and non-equilibrium 
electron transport through the double-island system 



shown in Fig. |l|. Structures consisting of two coupled 
metal islands or large quantum dots (large enough to 
blur the discreteness of energy JLavels) have recently re- 
ceived considerable attentiontHTtj. The present work dif- 
fers from the previous ones in two major respects. Firstly, 
in Refs. |ll|-|l4| and |l6| dots have been investigated which 
are coupled by just one or a few conducting channels, 
while in the present work a large number N of parallel 
channels in each junction is consideredO. (Ref. [15] also 
discusses how an arbitrary number of channels in the 
middle junction affects the positions of the conductance 
peaks.) Secondly, previous work has concentrated on two 
limiting cases: either strong tunneling between the two 
islands or between the islands and the leads, while we al- 
low for arbitrary relations between the various junction 
conductances. 
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FIG. 1. A schematic picture of a double-island structure 
with the voltage sources and various capacities in the system. 
The two circles denote the islands with ni /n r excess electrons. 

Systems of tunnel junctions are described by two char- 
acteristic quantities. The first is the charging energy for 
a given charge configuration (n\,n r ) on the islands, 
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It is displayed in Fig. |]for the four lowest states. The 
efficients -EVj&j with b = L, M, R, are the appropriate ca- 
pacitive energy scales, which in general depend on all ca- 
pacitances. The gate charges n^j = V g t,C g i,/e + VlCl/c 
and n x . x — V g B,C g Yi/e + VrCr/c control the charge on 
the islands, see Fig. || For convenience, we consider be- 
low a left-right symmetric situation: Cl = Cr = C, 

X)I = n x . The charging 
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energy contains two energy scales. One, 
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depends on the external capacitances and determines the 
curvature of the parabolas. The other scale, 
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also depends on the capacitive coupling, Cm, of the two 
islands and leads to a splitting as displayed in Fig. ^. The 
other important set of parameters are the dimensionless 
tunneling conductances through the junctions b 
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Also here we choose left- right symmetry, cvq = a^, 
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but 



do not fix the ratio between cvq 1 and a 




FIG. 2. Charging energy of the four lowest-lying states 
for the case of equal gate charges n x — n x \ = n XT . The states 
are denoted by (0,0), (1,0), (0,1), and (1,1). At temperatures 
T <C _Bcm/4, electron transport takes place in the regions 
around the degeneracy points of energy, depicted in the figure 
by the vertical dashed lines. 




FIG. 3. The honeycomb in the (n^i, n^^-plane shows the 
regions where a given state (ni, n r ) has the lowest energy. The 
circles denote the degeneracy points of Fig. ti. 



In this article we study the system using two ap- 
proaches. Firstly, we make use of the real-tinae diagram- 
matic technique developed for metal islandsBia and quan- 
tum dotsc3. This approach is applicable for linear and 
nonlinear response and all values of the gate voltages. In 
the main text, we discuss a systematic perturbative ex- 
pansion of the diagrams and, in appendix [A], present the 
calculation done up to second-order in the tunneling con- 
ductances. This calculation accounts for all second-order 
tunneling processes and turns out to be independent of 
any cutoffs. Secondly, in order to consider stronger tun- 
neling, we study the double-island system in terms .-of 
a recently-devoloped renormalization-group approachEj. 
In its general form, this procedure is cutoff independent 
and allows us to include all charge states. It thus differs 
from the RG approach developed in Ref. which only 
uses two adjacent charge states of the single-electron box 
and which requires a cutoff of the order of the charging 
energy Eq. We apply the RG equations in their equilib- 
rium form, which is sufficient for the linear conductance 
close to resonances, and consider tunneling conductances 
up to 8e 2 /h (a&'s up to 0.20). 

Our main results are the following. In the weak- 
tunneling limit, ajj < 1, and low temperature T <C 
Ecm/4: the linear conductance G displays a series of 
peaks as a function of the gate voltage. As the tun- 
neling conductances a b are increased, the peak heights 
are reduced with a temperature dependence resembling 
the functional form 1/[A + B\og(T/E c )} {A and B are 
T-independent parameters). Also the peak positions are 
renormalized due to higher-order processes. For the spe- 
cial case when the two islands are only weakly coupled 
to the leads, i.e., a^ 1 3> Q!q' R , we reproduce the a^f- 
dependences of the peak positions found in Refs. |l2| and 
rLq. As new features we find the corrections arising due 

— L R . 

to a ' including a weak temperature dependence of the 
peak positions. In the opposite limit, cvg' R ^-p^ 1 and 
just one channel in each junction, Andrei et al$3 predict 
a phase transition: the two conductance peaks around 
n x = 0.5 coalesce into one as Eqm is reduced below a 
finite critical value E^ M . By putting a^ 1 <C a^' R we ob- 
serve a similar trend for our junctions with large number 
of channels N. 

For completeness, we use the diagrammatic technique 
to account for the nonlinear conductance (dl/dV at V ^ 
0) and the inelastic cotunneling processes, which can- 
not be studied with the equilibrium RG procedure. The 
nonlinear conductance displays a double-peak structure 
resembling that found for SET's; in the present case the 
peaks are typically asymmetric. Second-order processes 
enhance this asymmetry and pull the two peaks further 
apart. On top of this, the inelastic cotunneling processes 
result in algebraically decaying tails of the peaks both in 
linear and nonlinear conductance. 
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II. MICROSCOPIC MODEL AND SEQUENTIAL 
TUNNELING 

The double-island system is described by the Hamilto- 



H = H +H T with 
Ho = Hl + Wr + Hi + H r + Tich- 



(3) 



The terms Hj describe noninteracting electrons in the 
leads and islands: Hj = ^2 k 



£ in a ]kn a 3kn with j £ 

{L, R} and Hi = J2 q , n e gn c iqn c iQn with i e {1, r}. The in- 
dices k and q label the electronic states, and n = 1 . . . N 
labels the transverse channels. We assume that the 
number of channels N in each junctions is large. The 
Coulomb interactions are described by the term 7Y c h 
which is the operator form of Eq. ([j]) . Tunneling of elec- 
trons is described in terms of the tunneling Hamiltonian 
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where T$ are the tunneling matrix elements for the bar- 
riers b (below we assume T^™ = T bn ). They are related 
to the tunneling resistances, e.g., for the right junction 
via 
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Here N R / T (0) is the density of states in the right 

lead/island and the operator e ±l< ^ changes the charge 
on the island i by ±e. 

Due to the periodicity of energies, it is sufficient to 
concentrate on gate voltages in the range < n x < 0.5. 
In this interval, the tunneling of electrons through the 
barrier b = L, M, R is characterized by the following dif- 
ferences in charging energies: Al = E c h(l, 0) — E cbj (0, 0), 
A R = E ch {0,l) - E ch (0,0), and A M = E ch (Q,l) - 
E ch (l,0). 

In the weak tunneling limit, cio < 1, sequential tunnel- 
ing processes dominante. At low temperatures only the 
states (0, 0), (1, 0), and (0, 1) participate in the transport 
and we obtain the following expression for the current 



jscq = 
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■-(A b ) = ±a b 



cxp [±/3(A b - A//6)] - 1 
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for the barriers b as well as the Fermi functions /^(A{,) = 
(exp[±/3(A(, — A/if,)] + l) _1 (A/if, is the difference in chem- 
ical potentials across barrier b). We further introduced 
a b = a+{A b )+a-{A b ). 

Already the sequential-tunneling current, Eq. (|6|), dis- 
plays interesting properties not found in single-island 
transistors. In the linear conductance regime one would 
expect to find the conductance peak at the degeneracy 
point, and this is indeed the case when all a\ are equal. 
In the asymmetric case, off ^ a ' , the peak acquires a 
temperature-dependent shift (vanishing at T = 0). The 
direction of the shift depends on the relative magnitude 
of aff and a$ = as shown in Fig. ^. Remarkably, 
the linear conductance at the degeneracy point is always 
given by the temperature-independent expression 
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corresponding to the intersection point, n x = 0.375, in 
Fig.|. 

In what follows, we account for charge fluctuations 
through all three junctions. We do this using two comple- 
mentary approaches: a real-time diagrammatic technique 
and a renormalization-group analysis. 




FIG. 4. 



Sequential-tunneling conductance for junc- 



tion 

«o = 
a^f = 0.01 



conductances a = a ' = 0.01 (solid line), 
= 0.001, ag' R = 0.01 (dotted line), and 
Oq' r = 0.001 (dashed line). The curves inter- 
sect at the degeneracy point of the charging energies due to 
scaling with Eq. (|^). For c^f larger (smaller) than a^' K in- 
creasing temperature shifts the peak to the left (right). 



For compactness of this expression, we have omitted the 
argument A^ from both the sequential-tunneling rates 
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III. DIAGRAMMATIC APPROACH 

The physical properties, including the nonequilibrium 
dynamics of the double-island system, are described by 
a reduced density matrix. We study the time evolution 
of this density matrix in a basis of charge states, \ — 
(n\,x> 7T. riX ), and perform an expansion in Ht. This leads 
to a formally exact master equation for the occupation 
probabilities p x in terms of general transition rates £ x , x ' 
which remain to be evaluated. In the stationary limit the 
master equation readalm 

= Px = bx' s x',x -Px s x,x'l ■ ( 9 ) 
x'#x 

Also the stationary current through the junction b can 
be expressed in terms of rates and probabilities, 

(I b ) = -ie^Px'^Xx- ( 10 ) 
x.x' 

The rate x with the plus ^ign is a subset of all pos- 
sible rates contained in S XiX 'Eil. We refer the reader to 
Refs. ^ and [l8] for a more thorough discussion on the 
above equation. 

The nonequilibrium time evolution of the density ma- 
trix may be visualized in terms of diagrams. In partic- 
ular, the rate ^ x , x > is a sum of all possible processes of 
which Fig. |B| shows some examples. In previous work, 
two approaches have been used for evaluating theutates. 
In the so-called resonant-tunneling approximationLm one 
occounts for an infinite sum of terms of a certain form 
via the use of self-consistent equations in the spirit of the 
Dyson equation. This approach accounts for high-order 
tunneling processes but is, in practice, limited to a small 
number of charge states. In the alternative, a systematic 
perturbative expansion of the diagrams^, one includes all 
processes of a given order O(a^). The key advantage of 
this approach is that, for not too high tunneling con- 
ductances, it can account for a larger number of charge 
states. 
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FIG. 5. Example of the diagrammatic representation of 
the time evolution of the reduced density matrix. In the di- 
agrams one can distinguish different tunneling processes (the 
vertical dotted lines are meant to guide the eye), from left 
to right, sequential tunneling, inelastic cotunneling, and sec- 
ond-order processes leading to vertex and propagator correc- 
tions. The pairs of numbers on the forward and backward 
propagators (the horizontal lines) correspond to charge states, 
the dashed lines with labels b describe tunneling of electrons 
across junction b. 



Thejarobabilities and rates can be expanded in powers 

S x,x' = E^, + S^, + . . . , ("J 

where p x and ^ xx i denote the terms ~ Uq. Figure |B| 

examplifies processes contributing to the rates T,^ with 
k = 1,2. When the expressions (|ll| ) are inserted into the 
master equation, Eq. (^), this must hold in each order 
and once we know the rates we may iteratively solve for 
p( k ) _ The probabilities and rates together determine the 
average charge on the islands (rii) — rii yX pi^ x as well 
as the current I = {lb), from Eq. (jlCj). 

A second-order calculation (O(oq)) that accounts for 
all charge states (reachable in this order) has proven its 
virtues in the study of SET's: the calculations could 
be carried out exactly and all cutoff dependences i£anish 
from the physical quantities. Also the comparisonQ with 
experiments showed that proper consideration of cotun- 
neling processes gives correct results for a wide range of 
tunneling conductances and also at resonances. A similar 
second-order calculation for the double-island system to- 
gether with rules for evaluating the individual diagrams 
is presented in appendix [A| Also the results for linear 
and nonlinear conductance are shown there. 



IV. RENORMALIZATION GROUP ANALYSIS 

In the strong tunneling regime, such as attained in 
recent experiments on single-electron transistorsa, it is 
necessary to go beyond low-order expansions in the tun- 
neling conductances. To this end, we study the double- 
island system using a recently developed real-time RG 
procedurelia. In its general form, this procedure accounts 
for an arbitrary number of charge states and is indepen- 
dent of any cutoffs. Below we use the equilibrium form 
of the RG equations which yields the renormalization of 
system parameters. We extend the equilibrium descrip- 
tion of the system to account for transport properties 
via the following assumption: for not too large the 
linear conductance is given up to small regular terms by 
the sequential-tunneling expressions when we replace all 
the bare parameters with the renormalized ones. The 
linear conductance obtained in this manner is consistent 
with the diagrammatic calculation at least up to O{o^). 
It should be noted that this approach does not account 
for the inelastic cotunneling processes. However, here we 
are mainly interested in the region around the conduc- 
tance maximum where, at least in the order O{o^), these 
only result in a weak temperature independent correction 
(similar result obtained for a SET). 

The results from the full RG are obtained numeri- 
cally and are thus presented graphically. For small a b 
the RG flow equations can be expanded analytically 
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up to 0(a^) and we can complement the numerical re- 
sults with the analytical formulas obtained in this or- 
der. The energy differences A^°\n x ) = -E c h(l, 0; n x ) — 
£ c h(0, 0; n x ) = E ch (0, 1 ; n x ) - E ch (0, 0; n x ) and the tun- 
neling conductances a\ are found to be renormalizcd 
to A(n x ) = A(°)(n x ) + A^{n x ) and d b = a b + a ( b 2) , 
respectively. The former yields directly the renormal- 
ized degeneracy point, n*, following from the condition 
A(n*) = A(°)(n*) + A^(n* x ) = 0, see Ref. |§ One 
should keep in mind that the peak arising from the ex- 
pression for the sequential-tunneling conductance is not 
necessarily located at the degeneracy point, see Fig. ||. 
The peak may be shifted due to unequal a\ but also for 
large Op because these can become significantly unequal 
in the renormalization flow. The full RG equations as 
well as their 0(aQ)-expansion (see below) for the general 
case are sketched in Appendix IB. 



V. RESULTS 

Figure ^ shows the main effects of increasing tunnel- 
ing strength on the linear conductance: the conductance 
peak is shifted and broadened and its maximum is re- 
duced. Furthermore the peak shape, in particular in 
the tail regions, becomes asymmetric. The key for un- 
derstanding the changes around the conductance peak 
maximum is the renormalization of the system param- 
eters: the peak lies at or close to the position of the 
degeneracy point, see Fig. 0, which is determined by 
the capacitive energies Eq and Eqm, the peak width 
depends on the curvature of the energy parabolas, Eq, 
and the peak height on the tunneling conductances db- 
For increasing tunneling strength, the effective capaci- 
tances C typically increase (modifying Ecb into Ecb) 
while db decrease. The fluctuations through the middle 
and the outer junctions affect the system in qualitatively 
different manner and below we distinguish between the 



three cases. 
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In all the figures, we have chosen left-right symmetry 
and, unless otherwise mentioned, used the convention 
Cl = Cr = Cm = lOOCgL = J.00CgR similar to the ex- 
perimental values of Joyez et aM. In terms of the charging 
energies, this corresponds to Ecm ~ ^cl,r ~ t-Ec- 
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FIG. 6. The solid line shows an example of the linear 
conductance peak as obtained with the full RG calculation. 
The dotted curve is the sequential-tunneling conductance for 
the same parameters. Here = Qo" = 0.025, a^f — 0.05 and 
\n(T/E c ) = -4. 

The case a^ 1 3> a L,R consists of a subsystem, the two 
islands, whose properties are probed via the weak tun- 
neling junctions connecting the subsystem to the leads. 
For weak coupling between the islands, i.e., small a^f, the 
electrons in each island form well-defined quantum states. 
For increasing a^ 1 , the states in different islands become 
increasingly mixed reducing the resulting ground-state 
energy. In terms of the system parameters, this can be 
seen as an increase in Cm - increase in Eqm - which low- 
ers the energy of the states (1,0) and (0,1) and shifts n* 
correspondingly, see Figs. and |^. The renormalization 
of cng's and A is independent of the gate voltage and the 
peak width remains given by thermal broadening. 
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FIG. 7. The position of the degeneracy point n% at T = 
as a function of Qq 1 with Oq' r = (dotted line, lower horizon- 
tal axis) and o^' R with ajf = (solid line, upper horizontal 
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Figure g displays the temperature dependence of the 
maximum conductance (the curves approach logarithmic 
form for off — ► 0). It shows also the maximum conduc- 
tance for c^f = 0.01 as obtained in second-order. The 
functional form is seen to differ from the correspond- 
ing RG result in three ways. First, the overall diverging 
shape is due to an unphysical /T-dependence of the 
perturbation expansion, emphasized for cases where the 
peak experiences a strong shift (large A^). Second, the 
initial slope at higher temperatures is slightly overesti- 
mated by the second-order result. Third, the maximum 
conductance is shifted downwards as compared to the 
RG result. This is due to the neglect of the inelastic 
cotunncling contributions in the RG calculation. 
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FIG. 9. The full conductance curve at ln(T/£ c ) = -5, -4 
for ag' R = 0.01 (dotted and dashed, respectively) and for 
ciq' r = 0.05 (solid and dot-dashed, respectively). The degen- 
eracy points of energy are located at 0.3385 and 0.3586 for 
c*o' R = 0.01, 0.05, respectively. 

Before going to the third case of equal ap's, let us write 
down the main results obtained by expanding the RG 
equations to O(o^) and make comparison to results in the 
literature. The renormalized degeneracy point, n*, fol- 
lows from the condition A«) = A^(n*) + A^(n* x ) = 
0. In the zero-temperature limit, the RG calculation 
yields the temperature- and n^-independent correction 
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for a n = 0.01, see text for further 



maxima for a^f- = 0.01, 0.05, 0.1. The curves approach log- For the special case E C l.r = E C m, the prefactors of a% 
0. The upward shooting curve is and a£ are equal. This linear ajf -dependence is a new 

result whereas the a^-dependence equals the results in 
Refs. |IJ,|l5|,|l2| (shown to be independent of N in Ref. jlB]) 
Further examination of the aj^-dependent curve in Fig. M 
in terms of polynomial fitting yields the coefficient of 
the quadratic term within 1-2% of the analytical result 
Ref. lM It should be noted that this accuracy already 
requires six correct decimals in calculating n*. The ojq- 
dependent curve in Fig. ffl remains very close to linear up 

L R 

to the largest a ' used here. However, if we repeat the 
calculation with a smaller value of Eqm the curve only 
approaches the line n x — 0.5 asymptotically. 

An expression for the temperature dependence of the 
peak conductance may be found by inserting the renor- 
malized a\ 's to Eq. (jg) . This corresponds to the require- 
ment that at resonance A(n*) = and yields 



arithmic form for a - 
the second-order result 
discussion on its shape. 

The limit a^' R 3> a^f corresponds to two inter- 
connected single-electron boxes (SEB) where the small 
conductance at the central junction can be used to probe 
the capacitive interactions of the two systems. Similar 
to single SEB's, increasing tunneling conductance in the 
junction b increases the capacitance Ch resulting in re- 
duced curvature of the energy parabolas and smaller A's. 
This can be seen as a broadening and shifting of the con- 
ductance peaks. The resulting peak shift is shown in 
Fig. [t] as a function of ccq ,r . Increasing tunneling conduc- 
tance also results in a gate-voltage and temperature de- 
pendent renormalization of (5l,r and 5m reducing these 
quantities. Figure ^ shows the full conductance curves 
for two choices of aQ ,R . The curious shape of the peaks 
is a consequence of the interplay of two effects: at a finite 
temperature the peak is shifted away from n* due to un- 
equal Qig's, see above, while ctb reduce the conductance 
around n*, i.e., not quite at the actual maximum. 
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with the temperature dependent parameters Bi. B± = 
2ln\T/E CM \, B 2 = 8\n\T/E CM \, and B 3 = ln\T/E CM \ 
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we obtain the limiting form 



G 



G^ cx -a 2 ln|T/^ CM | 
y 



(14) 



In this order, this is the sole source of temperature depen- 
dence to the maximum value of G. In the more general 
case cyq ' r ^ a^f , the peak is not at the point n* and one 
should use the full form for G. 

In the complete case, when charge fluctuations are al- 
lowed through all three junctions, the above-discussed 
effects intermingle in a complicated way. Now all the en- 
ergies Ecb and tunneling conductances a b are renormal- 
ized with some contributions to these depending on gate 
voltages and/or temperature while some others remain 
constant. In the following we consider for convenience 
the case of all a\ being equal. Even in this case, the 
peak may be shifted, cf. Eq. (|l^), but, in order to see 
the changes in the peak shapes more easily, the capaci- 
tances are chosen such that n* is independent of the tun- 
neling conductances. The choice of capacitances adopted 
above is close to the special case -EVjl.r = -E-cm with no 
shift in the second order and it appears that in this case 
there is no shift in higher orders either. Figure [lO] shows 
the full conductance curves for three values of a = 
at two temperatures in order to distinguish between the 
effects arising from thermal and quantum fluctuations. 
The broadening of the peaks with incresing ag is due to 
the overall reduction of A's as well as the gate- voltage 
dependent reduction of ctb- Figure ^ shows the tem- 
perature dependence of the maximum conductance for a 
few values of a®. At low temperatures, the curves resem- 
ble the functional form 1/(A + B m(T '/ E c )) , with the 
coefficient B being proportional to a$. This is consis- 
tent with the assumption that the full RG calculation 
resums the leading logarithmic terms ao\n(T/Ec)- Fig- 
ure [ll] also shows the second-order result for cvo = 0.01 
with the slope given by Eq. (|T|). The AW/T-behaviour 
of Fig. |s| is absent because the peak is not shifted and 
A««) = 0. 
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FIG. 10. The full conductance curve for ao = 

QfQ ,R =ao I = 0.01, 0.05, 0.10 (pairs of solid, dotted and dashed 
curves, respectively) calculated at ln(T'/i?c) = —4, —5 (upper 
and lower curve of each pair). The capacitances are chosen 
such that the peak does not shift either with temperature or 
with ao. 



G 
G 



1.1 



0.9 



0.7 



0.5 



0.2 




-9.0 



-7.0 



-5.0 



-3.0 -1.0 
ln(77£y 



FIG. 11. Temperature dependence of the maximum con- 
ductance for a = ao' R = off = 0.01, 0.01*, 0.03, 0.05, 0.10 
from up to down (* the curve with the diamond symbols is the 
second-order result for ao = 0.01). Below \n(T/Ec) ~ —3, 
the curves are of the functional form 1/(A + B\n(T/Ec))- 
The 



VI. CONCLUSIONS 

In conclusion, we have evaluated the transport prop- 
erties of a metal double-island structure using two ap- 
proaches: the real-time diagrammatic technique devel- 
oped in Ref. and the RG-procedure introduced in 
Rcf. |l9|. We made no a priori assumptions concerning 
the relative magnitudes of the tunneling conductances 
but we assumed a large number of channels in each junc- 
tion. The perturbative (diagrammatic) expansion has 
been performed up to O{o^), limiting the results to small 
values of ao. In this case we obtained both the linear 
and nonlinear response of the system for all values of 
gate voltage. In the RG calculation, we considered equi- 
librium properties of the system covering temperatures 
down to zero and also much larger tunnel conductances 
than in the perturbation theory. 

Due to periodicity of the charging energies, it is suffi- 
cient to concentrate on the single conductance peak ap- 
pearing in linear conductance in the gate-voltage range 
< n x < 0.5. The height and position of this peak are 
found to be renormalized and its shape becomes asym- 
metric. In general, we do not need to restrict the rela- 
tive magnitudes of the tunneling conductances ag, but 
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two limiting cases can be compared to existing work on 
similar systems. In the limit of chq 3> ocq — a§, we re- 
produce the third-order results Jor the peak position re- 
ported by Golden and HalperinEj. In the opposite limit, 
a^f <C «q = a^, the system resembles the model con- 
sidered by Andrei et alS3. These authors studied the 
mapping of a double-island system with one conducting 
channel in all junctions onto a Kondo model. It was pre- 
dicted that the system would undergo a phase transition 
as a function of Eqm , the mutual interaction energy for 
electrons in different islands: the two conductance peaks 
around n x — e/2 merge into one as Eqm is reduced be- 
low some finite critical value Eq M . By putting o^f small 
in our calculations for many-channel junctions, we inter- 
estingly obtain a similar trend of peaks shifting towards 
each other, see Fig. ^ 

We suggest the following experimental schemes similar 
to those used in Ref. || one could measure the heights 
both the linear and nonlinear conductance, as well as the 
asymmetric shape of the double-peak structure in the 
nonlinear conductace G(n x ,V 0). The tunneling con- 
ductance of a single junction is typically fixed in metallic 
structures. However, if one could attach an additional 
gate partially on top of a junction itself, one might be 
able to control the number of conducting channels instead 
of the conductance of separate channels. This would en- 
able the study of the peak position as a function of the 
tunneling conductances afj. 
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APPENDIX A: DIAGRAMMATIC 
CALCULATION TO 0(a%) 

In this appendix we present a perturbative calculation 
of the conductance through the double-island system up 
to the order O(o^). We also summarize the rules for eval- 

(2) 

uating the second-order rates ' , in the energy repre- 
sentation, i.e., in Fourier space. 

We begin by inserting both the probabilities and the 
rates in their series form, Eq. ( |TT| ) , into the master equa- 
tion Eq. (^J) and to Eq. ([it)]) for the current. We then pick 
up the terms of the orders O(ao) and 0{a^). The lowest- 
order contributions are of the form p^Y,^\ correspond- 



ing to sequential-tunneling processes, while in the second 
order two kinds of corrections appear. Terms such as 
p(o)5]( 2 ) correspond to two-electron tunneling -processes 
including the "inelastic cotunneling" processes! but also 
processes that give rise to renormalization of system pa- 
rameters. Terms of the formpWsW arise due to changes 
in the occupation probabilities in higher orders. While 
away from resonances the "inelastic cotunneling" yields 
the dominant contribution to the current, it was recently 
shown that close to the resonances it is crucial to consider 
all these second-order termsQ. 

The master equation, Eq. (^J) , must hold in each order 
and once we know the rates we may iteratively solve for 
p( k \ Hence, the essential step of the perturbative cal- 
culation consists of the evaluation of the transition rates 
X^ 1 ) and £( 2 ). The former are just the Golden rule rates 
and also the latter may be evaluated analytically as was 
done in Ref. |^. Before going into detail in evaluating 
the second-order rates, we would like to comment on two 
technical points. First, at intermediate stages we find it 
convenient to introduce a cutoff procedure to ascertain 
convergence of integrals. Any cutoff dependence cancels 
from all physical quantities, but only when we consider 
a set of 12 charge states - for the values of gate volt- 
age used here these are (0,0), (1,0), (0,1), (-1,0), 01-1), 
(1,1), (2,0), (0,2), (-1,1), (1,-1), (-1,2), and (2,-1)11 If 
we would not account for the additional charge states, 
the resulting quantities would contain terms with linear 
and logarithmic cutoff dependences. Second, for conve- 
nience, we solve the resulting set of 12 master equations 
numerically as this requires inverting a 12-by-12 matrix. 

1 . Rates 

The tunneling rates may be visualized in terms of dia- 
grams such as shown in Fig. [5], each diagram consisting of 
the following parts. The contour (solid line) corresponds 
to a forward and backward propagator, the upper and 
lower lines, respectively. Each part of the propagators 
has a charge state (ni,n r ) assigned to it and it carries 
the corresponding energy E nit n T . The directed dashed 
lines correspond to tunneling of single electrons through 
the junctions. These tunneling lines carry an index b, 
corresponding to the junction, as well as an energy lj of 
the tunneling electron. The diagrams, that cannot be di- 
vided into parts by a vertical cut through them without 
cutting one of the tunneling lines, are denoted irreducible 
diagrams and they correspond to rates of individual tun- 
neling processes. The tunneling rates between states x 
and \' (x — n hx' nr ,x) are obtained by summing up the 
rates of all the tunneling processes starting from the state 
X and ending up in the state x' ■ 

An arbitrary diagram may be evaluated as follows. 1. 
Each tunneling line corresponds to the expression 

a h (uj) — ±a n F —— rz -. (Al) 

b v ' °exp[±/3(w-A Mb )]-l v ' 
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with the plus (minus) sign corresponding to forward 
(backward) direction of the line with respect to the con- 
tour. The coefficient chq is the dimensionless conductance 
of the junction b, A/i& is the difference in chemical po- 
tentials across the junction and u> is the energy of the 
tunneling electron. 2. In order to account for all possible 
tunneling processes, one integrates over the energies u) 
and us' of the two tunneling lines. 3. The vertices, i.e., 
the end points of the tunneling lines, lying on the back- 
ward (lower) propagator acquire an extra factor —1 each. 
4. The diagrams can be divided into three parts i such 
that within each of them the energies of the propagators 
and tunneling lines remain constant. Each part i yields 
an energy denominator l/(A£'j + i0 + ) corresponding to 
the energy of the intermediate virtual state: AEi equals 
the difference of all energies of the lines (propagators and 
tunnclng lines) directed to the left and those directed to 
the right. The rate of a given process is the imaginary 
part of the expression obtained by evaluating the corre- 
sponding diagram. For example, we obtain for the third 
diagram in Fig. ^| 

- Mm • J dJa^) J dca+M ■ - - ^ + . Q+ 

1 1 (A2) 

uj + uj 1 + AE 2 +i0+ uj + AE 3 +i0+' v ' 

These integrals may be evaluated analytically along the 
lines of Ref. ||. It should be noted, that the use of these 
rules implicitly assumes a large number N of transverse 
channels in each tunnel junction. 

2. Results 

Conductance away from the peak is mainly given by 
the inelastic cotunncling processes and cannot be de- 
scribed by equilibrium calculations such as the RG anal- 
ysis in the main text. For certain values of the gate volt- 
ages, the lowest energy state is doubly degenerate, see 
Fig. [| In these regions we obtain a finite second-order 
contribution as, e.g., on the right-hand side of the peak 
in Fig. |l2|. The simplest processes that carry current 
through the whole double-island system are of the order 
0(o!q) and hence beyond the present considerations. For 
this reason we only obtain the exponential decay on the 
left side of the peak in Fig. [l2| Comparison between 
the present result and the RG calculation displays good 
agreement for small tunneling conductances thus sup- 
porting the assumption made (in RG) for the functional 
form of the conductance curve. 




-0.2 1 ■ ' ■ ' ^ 

0.36 0.41 0.46 „ 

FIG. 12. Sequential tunneling (dotted), cotunneling 
(dashed), and the total conductance (solid) as the sum of 
the two as obtained from the diagrammatic calculation. The 
logarithmic scale in the inset emphasizes the functional dif- 
ference of tails of the sequential-tunneling and total sec- 
ond-order conductances. The parameters used in the plots 
are Qq' r = af = 0.01 and T = 0.03E C . The junction ca- 
pacitances are now chosen unequal Cm = |Cl,r but still two 
orders of magnitude larger than the gate capacitances. This 
results in Ecm = |£7cl,r = f-Ec- 

A finite transport voltage applied across the system 
changes the charging energies as shown in Fig. At 
low temperature and when we only include sequential 
tunneling, the differential conductance dl/dV displays a 
double peak somewhat similar to that found in a SET 
but now the peaks are slightly asymmetric, see Fig. [TJ. 
For higher values of ao's the double-peak structure may 
be shifted as a whole, similar to the shift of the single 
peak in linear conductance. In order to inspect any ad- 
ditional effects resulting from the finite voltage, we have 
chosen the parameters in Fig. |l4| such that the center of 
the two peaks is not shifted. In this case, the peaks are 
shifted away from each other and the asymmetry of their 
heights and widths is enhanced. The different shapes of 
the tails are of the same origin as discussed above for the 
linear conductance, i.e., the asymmetric structure of the 
charging energies. 
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FIG. 13. Charging energy of the states (0,0), (1,0) and 
(0,1) when a finite symmetric bias voltage V is applied across 
the system, i.e., Vl = — Vr. The relevant current-carrying cy- 
cle of sequential-tunneling processes is depicted with arrows. 



^nun r = E ch{n\ + cr, n r ) - E ch (n h n r ), 

A « 1 CT « r = E ch {n h n T + a)- E ch {n h n r ), and (Bl) 

A^°n r = E ch(n h n T + cr) - Ech(ni + a, n r ) 

where a = ±1. The quantities corresponding to those 
used in the diagrammatic expansion are a^ x = a] 



■ 



a (b,o) = a o '> "(i,o) = a o > and A ( o,o) = Ab - 

We study the RG flow of the charging energies directly 
but the tunneling conductances are obtained from the 
auxiliary functions g defined as 



,Ti(0f-<M 



iL.na 



n L ± l,rc R )(ra L ,ra R | 
^ ^ T , a |n L ±l,n R Tl>(ni J ,n R | (B2) 



XI Sn^nR^L^R ± l)(n L ,n R | 
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FIG. 14. Nonlinear conductance through the dou- 
ble-island structure: sequential tunneling (dotted) and co- 
tunneling (dashed) contributions, and the total conductance 
(solid). The peak structure is asymmetric due to the asym- 
metry in charging energies with respect to the peak positions. 
The parameters in the figure are eV = 0.4_Ecm, T = 0.03-Ecm, 
E C m = £cl,r = %E C and all a£=0.01. 



APPENDIX B: RG FORMALISM 

In this Appendix we present central formulas of the 
RG approach applied irk-the second part of the text. The 
details may be found inl!3. 

We study the renormalization of charging energies 
£„ IjIlR and tunnel conductances a„ L „ R - As the sub- 
scripts already indicate, the various £J's and cv's renor- 
malizc in different ways for different charge states. For 
the following, let us define the differences of charging en- 
ergies more generally, 



These are introduced to parametrize the renormalization 
of «o's 



L 

HL,1R 

R 

iiL.nR 
M 

"L,™R 



a 9nt,n R 5n L + l,ri R 



a, 



R „R+ 



R- 



3HL,nR yn L ,riR + l 



(B3) 



a, 



M M + 



M- 



5n L ,n R 9n^ — 



Initially all g's equal unity corresponding to the bare a's. 

As an example of the full RG equations in the equilib- 
rium case we have for the energies 



dl llL.lR 



dl 



+ i 



^"•l,«r) — A n L ,n R 

i £(-l) x {<_ X+X ,,„ R e«p[(-l)x'iAj; L _ x+x , )nR ] 



"L+X,«R-X 

M 

nL+l-x+x'>nR-X 



_ Y ,ex P [-(-l)AA* 



"L + X.lR-X'J 



a 



_^exp[-(-l) x iA 



M 

«L + 1-X + X':«R-X 



(B4) 
} (B5) 



The summation indices x an< i x' take values zero and 
one. All quantities with a bar on top of them contain the 
time t c corresponding to the high-energy cutoff of the 
theory as, e.g., i£n L ,n R — E nijjnR t c . The parameter I is 
defined as I = ln(i c /i°) where t° c is the initial cutoff time. 
At the end, we take the limit t° c — > and integrate the 
equations to t c — > min{l/A, l/T}. 

The equations for charging energies are coupled with 
the flow equations for the „ R 's (cr = ±1) which follow 
equations such as 



d_ Lcr 

dl i<nL,riR 



E 

<t'=±1 



L Lcr' _Lcr L, — cr' 

^0 cr' ,nRynL-c' jJT'R^'MLjTT.r 



exp(iA 



i(A L<T l , - A L<T ' , ) 



10 



a 9ni J ,n R 9ni J -a' ,n R 9ni,n R ex P(*^n L £r',n R ) 

a R a Ra ' ,o Lct ,a R - a ' 

exp(iA^; g|nR _ g ,) - exp(zA%^ n _ ff ,) 
i(A RCT l , -A Rct ' ,) 

R L(j Rcr' R,— cr' „ fn A ^ ,< - T \ 

a 9n L ,n B .9ni J ,nB.—(7 , 9ni J ,n K ex P^^\i Lj n R — cr' J 

M Mcr' La /fl M,V 

ynL+CT+CT ,nR-j y^L+fT ,nR-(T "riL ,riR 

exp( ? A^; (T+ff ^ R _ (T ,) - exp(iA^; CT/! „ R _ CT ,) 



(B6) If A = the replacement 



z(A Mct 



nL+o-+er',nR— cr' 



A Ma ' , ,) 



a yn L ,n R 



o M<T ' , ,o M - <7 ' exv(iA Ma ' , ,)\ 



The full solution to these can be obtained via numerical 
integration. 

Expanding the above equations up to O{o^) yields the 
leading corrections to A's and a's. The resulting equa- 
tions are decoupled such that the equations for second- 
order quantities only depend on the respective bare quan- 
tities. For this reason the calculations can be carried out 
analytically. For the A's we arrive at 



_A fc '« 

dt r ™ L '" R 



(B7) 



— a e "l.»r - (n b -> n b+1 ). 

E ° b'cr 



The notation means that the second part of the equation 
equals the first part up to the indicated replacements of 
rVs. Here 6,6' £ {L, R} and A^ R = A%'_ a ,, nR - 
A^_ ff , . These equations can be solved to yield 

cxp( t n ;- nM 

, j A L<J (°) CiflA LCT(0) \ t) + A Lff(0) Si( A L<T(0) t ) 
-(n 6 ->n 6 + 1)1 (B8) 



The functions Ci(x) and Si(a;) are cosine and sine inte- 
grals. 

The functions g^(t c ) take the form as, e.g., 



> (1 \(tc)=$>0 



9n L ,n R ( 



'expHAl^)^)-^-^'^'^- 1 

A b'a'(0)~ 

LXfr, 



1 + aa' 



in L ,n R 



(B9) 



Ci (| A^'W + a*>c bb , | t c ) - Ci (| A^'W 
*Si ((A^'W + <7<t'c w ) i c ) + zSi ((A^'W) i, 

(*c-*2)]}. 



Ci^A^I^-T + ln 



Ab'a'(O) 



+ In(t c ) (BIO) 



has to be made. The summations over a and a' run over 
±1 and the parameters c bb i are given by 



cll = 2-Ecl, crr = 2-Ecr 
cmm = 2(i?cL + Eck — Ecm) 
cml = clm = £cm — 2i?CL 

clr = crl = Ecm 

cmr = CRM = 2I?cr — ^CM- 



(Bll) 



The quantities on the right-hand sides are the charging 
energies Ecb from Eq.(|l]). 
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